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Abstract 

Here, I first present the numerical computation of speed and direction of the drift of a spiral wave 
in an excitable medium in the presence of an electric field. In contrast to earlier results, the drift 
speed presents a strong variation close to the parameter value where the drift speed component 
along the field direction from parallel becomes anti-parallel. Using a simple phenomenological 
model and results from a numerical linear stability analysis of scroll waves, I show this behavior 
can be attributed to a resonance of the meander modes with the translation modes of the spiral 
wave. Extending this phenomenological model to scroll waves also clarifies the link between the 
drift and long wavelength instabilities of scroll waves. 
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Spiral waves can be observed in a variety of excitable systems such as Belousov- 
Zhabotinsky gels|li, colonies of the dictyostelium amoebae^ and slices of cardiac tissue^]. 
In the latter example, spiral waves of electrical activity have been shown to be the source of 
ventricular tachycardia and some of their instabilities are believed to be involved in the tran- 
sition from tachycardia to fibrillation, a deadly arrhythmia (For a review see ^y). This, with 
the intrinsic interest of those structures, has lead to an important research effort in order 
to understand the dynamic and instabilities of spiral waves and of their three dimensional 
analogous, scroll waves. 

In the presence of an external electric field, in the Belousov-Zhabotinsky reaction the 
center of rotation of spiral waves drift with a speed that presents components both parallel 
and perpendicular to the applied field The parallel component of the drift speed was 
always found to be in the direction of the applied field. A numerical study showed that 
depending on the parameter regime the drift direction of the spiral could be either parallel 
or anti-parallel to the field. The drift of a spiral wave has been linked 13, ^1 to the curvature 
instability of scroll waves I10| that leads scroll waves to bend and can finally result in a 
fibrillation like disordered activity of the medium. It was also linked to the three dimensional 
meander of scroll waves 11] which is the three dimensional analog of the meander instability 



characterized by a periodic modulation of the radius of rotation of the spiral wave. 
This phenomenon has been studied from a theoretical point of view ^, 0, . However, 

most analytical studies are restricted to the large core or the small core limit 1^ and 
can not examine the observed change in drift direction. 

In this paper, I present results of numerical computations that show an unexpected 
behavior of the drift speed at the drift direction change. The drift speed varies strongly in 
the vicinity of the transition. Using a reduced model of spiral wave in the presence of a 
small electric field, I show that this phenomenon can be attributed to a resonance between 
meander and translation modes of spiral waves. Finally, using the analogy between the 
effect of an electric field and the effects of a slight curvature of a scroll wave, I extend this 
model to scroll waves. This extension brings some light on the link between the drift of a 
spiral wave in an electric field and long wavelength instabilities of scroll waves and results 
obtained from this model are in good agreement with the results of the numerical linear 
stability analysis of scroll wavesQ]. 

I begin with describing the results of a numerical study of the spiral dynamics in the 
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FIG. 1: Tip trajectories in the presence of an electric field for a equal (left to right) to 0.85, 0.91, 
0.93 and 0.99 for E=0.03 (top) and E=0.003 bottom (Other parameter values are h = 0.13 and 
e = 0.02). The trajectories were recorded for 100. time units (tu) (top) and 400 tu (bottom) 
and are represented in boxes of same spatial extension. The arrows show the drift direction. For 
parameter values a = 0.85 and a = 0.99, the drift velocity behaves almost linearly with the strength 
of the field (this was checked by performing two simulations, see also fig. |2l (one should note that 
this is enough to check wether the linear regime is reached since an additional point is f = 0, 
E = 0) while in the cases a = 0.91 and a = 0.93 the behavior is still strongly non-linear. 



presence of an electric field using tlie Barkley model of an excitable medium (contrary to 
BZ systems, u is the only diffusive variable): 

dtu = -f{u,v) + Au + EdxU, (1) 
dtv = g{u,v), (2) 

where f{u, v) = u(l — u){u— ^^), g{u^ v) = u — v and E is the electric field, directed along 
tfie X axis. An homogeneous system modeled by these equations presents a single stable 
equilibrium point u = v = Q and a small perturbation of this equilibrium point can lead to 
a large excursion in phase space before return to equilibrium. The diffusive coupling in Eq. 
fallows the propagation of solitary waves in ID. In 2D, in the absence of an electric field 
(-E = 0), a rich variety of wave propagation regime is observed depending on the values of 
the parameters jl?! . I will focus here on steadily rotating s pira ls and present results obtained 
along a line of equation 6 = 0.13 in the parameter space of 
spiral is observed). 
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In the presence of the electric field {E ^ 0), the spiral tip drifts with a constant veloc- 
ity (see figure that depends on the parameter values and on the electric field strength. 
The drift velocity has components perpendicular {v±) and parallel (fy) to the field. These 
velocities vary linearly with the field in the small E limit. This allows to define the drift 
coefficient a± = v±/E and a\\ = v\\/E^^ in a weak field. 

Numerical results (see Fig. ^ show that close to the transition from anti-parallel to 
parallel drift, the value of E for which the linear regime is reached is much smaller than for 
the small a and large a region. In addition, close to this transition, the dependence of the 
drift coefficients a\\ + ia± in a is strongly non monotonous. When increasing a there is a 
strong enhancement of the parallel drift coefficient followed by a relatively sharp transition 
from parallel to anti-parallel drift and finally the drift coefficients decrease rapidly. The 
behavior of the perpendicular drift coefficient is characterized by a strong enhancement 
close to the value of a where the transition occurs. These results differ strongly from the 
results presented in where the drift speeds were computed using the same parameter 
values and an electric field of amplitude E = 0.3: a monotonic variation of the drift speeds 
with the parameter values which is also observed here when using high values of E for which 
(see Fig. |21 solid line) the drift amplitude is no longer linear with the field amplitude. 

I now present an ODE model of spiral wave drift in an electric field that explains the 
strong enhancement of drift coefficients by the resonance of damped meander modes and 
translation modes. 

The ODE model is a modification of models 0, Q] that reproduces the main feature of 
the spiral wave dynamics and has been modified in order to take into account the effects of 
a small electric field 

f = e'^iRoUt -z)+(3E (3) 
z = [^i-iujm-{l+ ia)\z'^\]z + ^Ee-'^ (4) 
(f) = uJt (5) 

where r is the position of the tip in the complex plane ((x, y) in the real plane correspond to 
(i?e, Im) in the complex plane) and 2; is a complex variable describing the meander in the 
frame rotating with the spiral and yU + ium is the eigenvalue associated with the meander 
mode and E is the amplitude of the field directed along the real axis in the imaginary plane. 
In the absence of an electric field (i? = 0), for < 0, the steady state is obtained for z = 
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FIG. 2: Parameter values are e = 0.02 and b = 0.13. Top left: v^^/E for the following values 
of the applied external field: 0.3( solid), 0.1 (long dashed), 0.01 (short dashed), 0.001 (dotted), 
0.0001(+), the dash-dotted line is the line of equation v^^/E = 0. Top right: v±/E for the same 
values of the external field, E. Note that for both larger and smaller values of a (not shown here) 
v^^/E and v±/E are independent of E for the E range explored here. Note that the curves for 
E = 0.0001 and E = 0.001 are indistinguishable, indicating that for those values, the drift speed is 
linear in E. Hence the dotted curve also represent a± and ay. One should note the very different 
shape of the two dotted curves. Similar results (with a not so strong drift enhancement) were 
observed for higher values of b and stronger damping of meander modes (6 = 0.18). Bottom left: 
schematic of the phase diagram in Barkley's parameter space. In the region left from the dotted 
line no spiral waves are observed. On the right hand side of this line, rotating and meandering 
spiral waves (gray region) are observed. The dashed line corresponds roughly to cOm = and the 
horizontal solid lines correspond to the values of b (0.13 and 0.18) for which systematic simulations 
were performed. Bottom right: as a function of log^g ^ — 0-92 and b = 0.13. E is varying 
from 0.0001 to 0.3 and the dashed line is of equation t!|| = 9.3t;. 



and corresponds to the stationary rotating spiral r = i?o exp (icUft) at frequency Ut. For 
yU > 0, the z = solution of Eq. E]is no longer stable and the z variable undergoes a Hopf 
bifurcation that leads to the modulation of the radius of curvature of the tip with frequency 
—ujm characteristic of meandering spirals (one should note that for Ut = ujm, the meander 
instability (/i > 0)leads to the drift of the rotating spiral with constant velocity [17] for 
E = 0). 
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I now give some rationale for the terms introduced in eqs. ()3l4l5p . In eq. El the effect of 
the field can be either independent of the relative orientation of the spiral and the field (the 
f3E term) or depend on it. Nonetheless, the later case leads to a precessing term that once 
integrated over one spiral rotation vanishes at dominant order. Hence, I have chosen to add 
only the former term and to disregard the later that would have no effect at leading order. 
The 'f E exp{—i(j)) term in eq. I^has to take into account that the field is constant in the 
laboratory frame while eq. |3] is expressed in the frame rotating with the spiral, that is the 
laboratory frame rotated of 0. Therefore it is necessary to add the exp(— i0) factor which 
expresses the fact that the meander mode is affected by the field in a way that depends on 
the relative orientations of the spiral and of the electrical field (in the frame rotating with 
the spiral the field is rotating with a pulsation which is —0). No extra term has been added 
to Eq. El (also expressed in the frame rotating with the spiral) since at dominant order the 
addition of a a term proportional to E exp{—i(j)) would only lead to another constant drift 
term in Eq. 121 similar to the f3E term. 

I focus here on the case where steadily rotating spirals are stable, (i.e. /i < and ov is 
the frequency of the meander mode that can be computed by linear stability analysis j20| of 
steady spirals) and describe the effects of a small field in this situation. At the leading order 
in E and z, Eq. HI has for solution z = Aexp{—iutt) with A = 'yE/{—fi + i{—ujt + ojm))- 
Using this expression in Eq. El leads to a spiral tip drifting with constant velocity: 

^ ^ ^ , 7(/^ - i{uJt - UJm)) 

E f,^ + (uj, - um)' ^ ^ 

This expression describes qualitatively the behavior of the spiral drift coefficients in the 
presence of an electric field as a function of a presented in Fig. |21 if the two following 
conditions are met. First close to the point where the transition from anti-parallel to parallel 
drift occurs the meander frequency, ujm becomes equal to the spiral frequency Ut and fi is 
small {weakly damped meander regime). Second, the real part of 7 is much smaller than its 
imaginary part in this region. Else, the expression of Eq. El can not reproduce qualitatively 
the behavior observed in fig. El One can check this statement looking at the value of the 
drift coefficient in the case where 7 = i'-^i is imaginary. In this case, the expression of the 
drift speed is 

E jj,^ + {uJt - iO^y /i2 + {uJt - 

which reproduces well the behavior presented in fig. El when n is small. 



It is also interesting to note that the drift coefficient as expressed in eq. p diverges if 
/i = and Um = ^t, that is at the codimension two point of Barkley phase space |l7| (noted P 
here). This seems unphysical since in the absence of an electric field one can see a perfectly 
stable spiral. Nonetheless, as shown in for Um = and /i > the meandering spiral 
tip trajectory is the one of a steadily drifting spiral with a constant finite speed even if 
E = Q (this situation corresponds to the case of an infinite drift coefficient). Hence it is not 
surprising that the drift coefficient should diverge when approaching P in the stable spiral 
region. 

The following part of this paper will show that these requirements are met. In addition, 
the link between the drift of a spiral wave and three dimensional instabilities of scroll waves 
will be discussed. 

One can determine the values of Ut and Um using the linear stability analysis of steady 
spiral previously described in [2^ and extended to scroll waves in As seen in Fig. El oJt 
and ujm become equal for a 0.925 (i.e. close to the transition from anti-parallel to parallel 
drift). In addition, for this value of a, the damping of meander mode is relatively weak. 
Finally, taking 7 = 0.4^ and (3 = —0.36 — 1.3i, the behavior of the drift coefficients is well 
approximated by the expression of Eq. ^ in the vicinity of the transition. 

Since the effects of a small electric field on a spiral are analogous to the effects of the 
slight curvature of a scroll wave, it is of interest to consider the three dimensional extension 
of this model. In the following part of the paper, I will describe the results obtained using 
this simple model. The results presented here conffim that the low kz curvature of the 
translation modej^ is equal to the drift coefficients. They will also show that contrarily to 



what was hypothesized in the drift coefficients are not equal to the opposite of the low 
kz curvature of meander modes. 

First, consider a slightly curved spiral ffiament. The evolution equations of u and v are 
the three dimensional analog of Eqs. ^121 in the case E = 0. They can be rewritten using 
the coordinates s, X, ,Y where s is the curvilinear abscissa along the spiral ffiament 21|, 
X and Y are the cartesian coordinates in the plane perpendicular to the ffiament (the X 
axis being oriented along the normal to the filament). Using this coordinate system and 
assuming a slight curvature of the filament, Eq. |21is unchanged and Eq. breads: 

dtu = -f{u, v) + Axyu ^-—^dxu, (8) 

6 p — A 
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FIG. 3: Top left: computed frequencies of the meander mode (x) and of the translation mode 
(+) as a function of the control parameter a (e = 0.02 and b = 0.13). uj^ and ujt are equal for 
a ~ 0.925. Top right: growth rate of the meander mode as a function of a (other parameters 
unchanged). The maximal value of /u is reached for a ~ 0.920. Bottom left: curvatures of the 
real part of the translation branches (+) and opposite of the curvature of the real part of the 
meander branches (x) at /c^ = 0. Bottom right: Computed curvatures of the imaginary part of the 
translation branches and opposite of the curvature of the imaginary part of the meander branches 
at kz = 0. For a = 0.92 and 0.93 the computed values of the curvatures present an uncertainty of 
order 1 due to the hybridization of meander and translation modes. Dashed line in bottom figures: 
fit obtained using Eq. ^Jfor the translation modes and a second order polynomial fit of 7(a) and 
P{a). Fit for the meander modes not shown here. 



where p is the radius of curvature of the filament and Axy is the two dimensional laplacian 
operator in the XY plane. The — (l/(p — X j)dxu term, at leading order in X/p is equal to 
1/p. Hence in the new frame, close to the filament, Eq. ^is affected by the curvature of the 
filament in the same way it is affected by an electric field E = —1/p. 

Therefore, the three dimensional extension of the Barkley-Sandstede model for slightly 
curved filaments is given by eq. ()3I4I5|) where E is replaced by — r" and where r" denotes the 
second order derivative of r along the axis of the filament. For slightly curved filaments r" is 
equal to the opposite of its curvature. A linear stability analysis of the tridimensional ODE 
model around the steady scroll wave {z = 0, ro(-2,t) = exp{iujtt)) in the frame rotating with 
the spiral i.e. using a perturbation of the form r = ro(-2, t) + 7--^e«0e'^*+«'=zz g^^^ ^ _ ^_^^crt+ik^z 
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leads to the following eigenvalue problem: 



Zl 




(9) 



whose eigenvalues are at leading order in k^: 



a, = -tuj, + + kl (10) 

(T^ = /i - iu^ - I — —P- ) kl (11) 

These equations show that the curvature of the translation branch is equal to the opposite 
of the curvature of the meander mode shifted of (3. It is also equal to the drift speed (see 
Eqini). Hence, contrary to what was proposed by Aranson et al. 11|, the curvature of the 
meander mode is not equal to the drift coefficients. There is a shift (that can be relatively 
small when the parameter values are close o the line where ujm = cuj as shown in Fig. 01 
This might come from the fact that the parameter P was omitted in jll']. 

In addition, since the numerical linear stability analysis of scroll waves allows to determine 
at, <Jm, /U, Ut, and Um for small values of kz, one can extract from Eqs ilOllH the expressions 
of the coefficients /5 and 7: 

/3 = -{k^ + Kt) (12) 

7 = KmijJ^ + i{uJt - i^m)) (13) 

• Urn) + {Vm • Vm)) (14) 
Kt = ~{Ut»Ut)/i{Ut»Ut) + {Vt»Vt)) (15) 

where Km and Kt are the respective curvatures of meander and translation branches at k^ = 
and can be expressed as functions of {um,Vm), {ut,Vt) the meander and translation modes. 



{um,Vm), {ut,Vt) are the corresponding adjoint modes that can be computedj22|| and (•) 
denotes the usual scalar product in the frame rotating with the spiral. 

This method in the vicinity of uJm = oJt results in great variations in the values of 7 since it 
is very sensitive to possible inaccuracies. However, for values of a away from the transition, 
one can compute with a good accuracy both /5(a) and 7(a). They appear then to be smooth 
functions of a that can be easily fitted by a second order in a polynomial. For a close to 
the transition from parallel to anti-parallel drift, the polynomial fits of (3{a) and 7(a) takes 
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values close to the ones used previously to fit the drift speed and using the polynomial fit 
in Eqs. ITUllTTl one reproduces with a good accuracy the results of the long wavelength linear 
stability analysis of scroll waves (see Fig. |2)). One should also note that Eq. ITHl together 
with the results of the numerical three dimensional linear stability analysis od scroll waves 
give some rationale for 7 being imaginary close to the transition from anti-parallel to parallel 
drift. Indeed, close to the transition Km is mainly directed along the imaginary axis in the 
complex plane and Um ~ ^^t, those two results lead to 7 being directed almost along the 
imaginary axis. 

To conclude, using smoothly varying coefficients (3 and 7, this model quantitatively re- 
produces the results obtained numerically when considering both the drift of a spiral wave 
in the presence of an external field and the long wavelength instabilities of scroll waves (me- 
ander and curvature). The reduced model brings some clarification on the mechanism of 
the drift of a spiral wave in the presence of an electric field, showing that the change in the 
sign of drift velocity parallel to the electric field can be attributed to a resonance between 
meander and translation modes and that this resonance leads to a strong increase in the 
drift coefficients. One should also note that despite the results presented here are formally 
valid for /i << 1, the resonance described here infiuences the drift speed as a function of 
parameters for a wider parameter regime and its effects can be observed rather far away 
from the Barkley's codimension 2 point (see fig. 2 a (inset)) as long as meander modes 
are not too damped. I also hope those results can be checked exp erimentallv in BZ system 
where meandering of spiral waves has been well characterized |23l| . 
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